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Abstract 

We study the spin Calogero model oi Dm type with polarized spin reversal operators, 
as well as its associated spin chain of Haldane-Shastry type, both in the antiferromag¬ 
netic and ferromagnetic cases. We compute the spectrum and the partition function of 
the former model in closed form, from which we derive an exact formula for the chain’s 
partition function in terms of products of partition functions of Polychronakos-Frahm 
spin chains of type A. Using a recursion relation for the latter partition functions that 
we derive in the paper, we are able to numerically evaluate the partition function, and 
thus the spectrum, of the Z3jv-type spin chain for relatively high values of the number 
of spins N. We analyze several global properties of the chain’s spectrum, such as the 
asymptotic level density, the distribution of consecutive spacings of the unfolded spec¬ 
trum, and the average degeneracy. In particular, our results suggest that this chain is 
invariant under a suitable Yangian group, and that its spectrum coincides with that of 
a Yangian-invariant vertex model with linear energy function and dispersion relation. 

Keywords: Exactly solvable quantum models, Spin chains with long-range 
interactions. Polarized spin reversal operator 
PACS: 02.30.Ik, 75.10.Pq, 75.10.Jm, 05.30.-d 


1. Introduction 

In a recent paper, a novel type of spin Calogero models and their associated spin 
chains of Haldane-Shastry type was introduced [T]. The distinguishing feature of these 
models is that they are constructed using a new representation of the Weyl group of 
the BCn root system, obtained by replacing the standard spin reversal operators by an 
arbitrarily polarized version thereof. As shown in the latter reference, these models are 
exactly solvable for all such representations and, in particular, the partition function of 
the spin chains can be exactly computed using Polychronakos’s freezing trick W- In this 
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paper, we shall extend the above results to spin Calogero models and their corresponding 
spin chains of Haldane-Shastry type based on the Djq root system. 

In order to present our work in an appropriate context, let us briefly recall the origin 
and significance of the latter models. The Haldane-Shastry (HS) spin chain, introduced 
independently by these authors in the late eighties mis], is perhaps the best known ex¬ 
ample of an exactly solvable one-dimensional lattice model with long-range interactions. 
More precisely, this model describes a circular array of equispaced spins with two-body in¬ 
teractions inversely proportional to the square of the (chord) distance between the spins. 
The motivation for introducing this chain was the construction of a simple model with 
an exact ground state given by the C/ —>■ oo limit of Gutzwiller’s variational wavefunction 
for the ground state of the one-dimensional Hubbard model [SHE]. Over the years, the 
HS chain has appeared in many areas of interest both in Physics and Mathematics, such 
as fractional statistics and one-dimensional anyons inHH], quantum entanglement m, 
characterization of integrability vs. quantum chaos HMD, quantum integrability via the 
asymptotic Bethe ansatz [THHEO], Yangian quantum groups HaillHESl, and conformal 
field theory [181 124H26] . 

One of the key properties of the HS chain —already noted by Haldane and Shastry in 
their original papers— is its intimate connection with the scalar (trigonometric) Suther¬ 
land model HZlllH]. This connection was subsequently elucidated by Polychronakos in 
Ref. |2] , who showed how to derive the HS chain from the spin Sutherland model 
by a technique that he called the “freezing trick”. The main idea behind this technique 
is to note that when the coupling constant in the spin Sutherland model goes to infinity 
the particles tend to concentrate on the coordinates of the (essentially unique) minimum 
of the scalar part of the interaction potential, which are precisely the sites of the HS 
chain. Thus, in this limit the dynamical and the spin degrees of freedom decouple, and 
the latter are governed by the chain’s Hamiltonian. Using this idea it is straightforward, 
for instance, to obtain the first integrals of the HS chain from their well-known counter¬ 
parts for the spin Sutherland model. In fact, Polychronakos showed that applying the 
same procedure to the (rational) Calogero model [32] and its spin version |3T] one ob¬ 
tains an integrable chain with non-equispaced sites and long-range interactions inversely 
proportional to the distance between the spins |2] . The spectrum of this chain —known 
in the literature as the Polychronakos-Frahm (PF) chain— was first studied numerically 
by Frahm [33j and then exactly computed by Polychronakos |3], who derived an exact 
formula for the partition function by means of the freezing trick. On the other hand, the 
partition function of the HS chain was only evaluated more than a decade later by some 
of the authors HH. 

Both the Sutherland and Calogero models (and their corresponding HS and PF spin 
chains) are associated with the Ajv-i root system, where N is the number of particles. 
Indeed, in these models the interactions only depend on the difference between the coor¬ 
dinates, and the spin operators appearing in the Hamiltonian are permutation operators, 
and thus generate a realization of the Weyl group of type. In fact, there are ver¬ 

sions of the Sutherland and Calogero models associated to any (extended) classical root 
system |34|. Among these systems, those of BCn, Bn, Cn and Dn type are by far the 
most studied in the literature, since they make it possible to construct integrable models 
with an arbitrary number of particles. By applying the freezing trick to the spin version 
of these models one obtains the corresponding generalizations of the HS and PF chains, 
that we shall collectively refer to as spin chains of HS type [331441] . One of the funda- 
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mental features of the BCj^ root system and its Bjv and Cjv degenerations is the fact 
that its Weyl algebra contains a family of reflection operators Si {i = 1,..., N) satisfying 
Sf = 1. (In the case of the root system, the Weyl group only contains products SiSj 
with i ^ j.) In the spin chains studied in Refs. [3SJH], the operators St are represented 
by spin reversal operators Pi (acting on the Hilbert space of the i-th particle), but this 
is by no means the only possible choice. As a matter of fact, in the novel version of 
the spin Calogero model of BCn type and its corresponding (PF) chain introduced in 
Ref. [T], the operators Si are represented instead by arbitrarily polarized spin reversal 
operators (PSRO) which act as the identity on the first mi elements of the 

spin basis and as minus the identity on the rest. These operators are equivalent under 
a similarity transformation to the usual spin reversal operators Pi only when mi = m 2 
or mi = m 2 ± 1, i.e., when there is minimal polarization. For the remaining values 
of the discrete parameters mi and m 2 , the systems constructed in the latter reference 
differ from the standard Calogero and PF models of BCAr-type. In particular, when mi 
or m 2 are zero, the corresponding spin dynamical model reduces to the su(m)-invariant 
extension of the Calogero model studied by Simons and Altshuler (42] ; see also [43] . 

In this paper we introduce the spin Calogero model of Dp^-iy^e with PSRO and 
its corresponding spin chain of HS type, i.e., the PF chain of Dp; type with PSRO. 
As explained in Ref. |3S|, these models are singular limits of their corresponding BCn 
counterparts, so that their spectrum cannot be obtained by setting to zero the parameter 
P in the latter models (cf. Eqs. (2.5) and (2.14)). This is also apparent at the level of 
the Hilbert space, which is the direct sum of the Hilbert spaces of two BCn models with 
opposite chiralities. Thus, the models studied in this paper are not limiting cases of their 
BCn versions in Ref. [T]. 

Our main result is the derivation of a closed-form expression for the partition function 
of the PF chain of Dn type with PSRO in terms of products of partition functions of 
type-A PF chains. Our approach is based on the computation of the spectrum and 
partition function of the corresponding spin Calogero model, from which the chain’s 
partition function follows by a standard freezing trick argument. The structure of this 
partition function turns out to be more involved than that of its BCn counterpart. In 
particular, it is not manifest that it is a polynomial in g = as follows from 

the freezing trick. Using the explicit expression for the partition function, we shall study 
several global properties of the chain’s spectrum, such as the behavior of the level density 
and the average degeneracy when the number of spins tends to infinity. In particular, 
the fact that the number of distinct levels grows polynomially with the number of spins 
suggests that this model is isospectral to a Yangian-invariant vertex model of the kind 
studied in Ref. jUj. 

The paper is organized as follows. In Section we recall the definition and main 
properties of the polarized spin reversal operators construct the Hamil¬ 

tonians of the Ujv-type spin Calogero model with PSRO and its associated spin chain. 
Section [^ is devoted to the derivation of the closed-form expression of the chain’s parti¬ 
tion function, as explained in the previous paragraph. Using this expression, in Section]^ 
we analyze several global properties of the spectrum, providing strong numerical evi¬ 
dence of the Gaussian character of the level density when the number of spins is large 
enough. In Section we extend the above results to the ferromagnetic version of the 
models under consideration. The paper ends with a brief section summarizing our main 
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results and presenting our conclusions, and a short technical Appendix establishing a 
useful recursion relation for the partition function of the PF chain of type Ajv-i- 


2. Construction of the models 


For the purpose of describing the Z?Ar-type Calogero model with polarized spin reversal 
operators, it is convenient to briefly summarize the construction of its BCjq counterpart 
[T]. To this end, let 

5 = (|si,... ,SAr) |si G {l,2,...,m}). (2.1) 

denote the internal spin space for N particles. As usual, the action of the spin exchange 
operator Pij on S is defined as 


••• T 1 -i ) — I • 


, CJj , . . . , 02 


,..., sn) ■ 


( 2 . 2 ) 


Let us denote the PSRO associated with the Fth particle as where mi and 

TO 2 are two nonnegative integers satisfying the relation mi + m 2 = m. The action of 




where 


on S is given by [I] 


p(mi,m2)| 

-Ci *1) 


f{s^) = 


,■■■, Sn) = (—l)'^^^’^|si5.. •, Si,..., Sn) 

f 0, 1 ^ Sj < TOi 

I 1, mi + 1 ^ Sj < TOi + m 2 . 


(2.3) 


(2.4) 


In terms of these operators, the Hamiltonian of the HCiv-type Calogero model with 
PSRO is defined as [T] 


H, 


{mi,m2) 
B,e 






Ct + Pij 


p(mi,m2) ' 


(4-)^ 


+/3a^ 


/3a - eP- 


{mi,m2) 


(2.5) 

where the sums run from 1 to A^, a > |, /3 > 0, e = ±1, a:)':- = Xi± Xj, = X^i 


P 


(mi,m 2 ) 


_ p{mi,m2) p{mi,m2) p 


ij * 1 1 * *1 ■ 

It can be shown that when m is even (resp. odd) and mi = m 2 (resp. mi = m 2 ± 1), 
the PSRO in (2.3) is equivalent via a similarity transformation to the usual spin reversal 
operator Pi, which changes Si into m — Si + 1. As a result, for these special choices of 
mi and m 2 , the Hamiltonian (2.5) reduces to that of the standard su{m) spin Calogero 
model of BCn type studied in Ref. [35] • As mentioned in the Introduction, another 
interesting special case is mi = m, m 2 = 0, for which the Hamiltonian ( |2.5| ) reduces to 
the Simons-Altshuler extension of the Calogero model. 

Since contains the discrete parameters mi, m 2 and e, it is natural to inquire 

whether there exists any relation between models (2.5) with different sets of parameters. 
In fact, we shall now show that is equivalent to through a unitary 

transformation. To this end, consider the unitary operator T whose action on the spin 
space S is given by 

P|Sl,...,Sl,...,S7v) I'^l5---5'^2J---J S N ) ; 
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( 2 . 6 ) 








where 


Si + mi, 1 ^ Si ^ m 2 

Si — m 2 , TO2 + 1 < Si < mi + m2- 


Using Eqs. ( |2.3[ ) and (2.6) we easily obtain 

T^Pf"")T\si, ...,s„...,sn) = (-l)^(*-^|si, ...,s,,...,sn). 


(2.7) 


( 2 . 8 ) 


From Eq. (2.7) it follows that s' G {mi + 1,..., mi + m 2 } for 1 ^ Si < m 2 and s' G 
mi} for m 2 + 1 ^ Si ^ m 2 + mi, so that 

/(s') = I 


1 < Sj < m 2 

m2 + 1 ^ Si < m2 + mi. 


(2.9) 


Equations (2.8) and (2.9) clearly imply that 


It is also obvious from Eqs. (2.2) and (2.6) that 

P T=P 

From Eqs. (2.5), (2.10) and ( |2.11[ ) we readily obtain 

pj-(mi,m2)p _ jp{m2,mi) 


( 2 . 10 ) 


( 2 . 11 ) 


( 2 . 12 ) 


as claimed. In view of the above relation, it suffices to study the Hamiltonian (2.5) in 
the case e = 1. However, in the paper we shall intentionally keep the parameter e in 
in order to facilitate the comparison with its Ujv counterpart that we shall 
introduce below. 

Due to the nature of the singularities of the Hamiltonian configuration 

space can be taken as one of the Weyl chambers of the BCm root system, i.e., one of the 
maximal open subsets of on which the functions Xi ± Xj and Xi have constants signs. 
We shall choose this configuration space as the principal Weyl chamber 


= |x G : 0 < xi < X 2 < ■ ■ ■ < a;Ar} , 


(2.13) 


where x = {xi,... ,xn). The Hamiltonian jg thus defined on an appropriate 

dense subset of the Hilbert space ® S. When e = 1, the spectrum of 

was computed in Ref. [T] by constructing a suitable (non-orthogonal) basis of this Hilbert 

space in which this Hamiltonian acts triangularly. 

As explained in the latter reference, from the spin dynamical model (2.5) one can 
construct a PE chain of BC^ type with PSRO by applying the freezing trick. The 
Hamiltonian of this chain is given by 


n 


{mi,m2) 


E 


1 + Pij 

(G-0)^ 


i + Pi 


(mi,m2) 


(Ci + 0)^ 


'E 


1 - eP 


(mi,m 2 ) 


Si 


(2.14) 



















where the lattice sites Q are related to the zeros of the Laguerre polynomial ^ by 
yi = Q/2. The exact partition function of the chain (2.14) has also been computed in 


Ref. [T] by exploiting its connection with the spin dynamical model (2.5). 

The Hamiltonian Qf Ujv-type spin Calogero model with PSRO is nat¬ 
urally defined by dropping the term related to the roots Xi in ^ j^g,^ by setting 

/3 = 0 in Eq. (2.5). We thus obtain 




P,. 


a + py 


mi ,m2) 




(4)^ 


(2.15) 


It should be noted that, unlike its BCff count erpart, the latter Hamiltonian does not 
depend on e. Just as before, from Eqs. (2.10) and (2.11) it follows that jg 

unitarily equivalent to under T: 




(2.16) 


Thus we can impose without loss of generality the restriction mi ^ m 2 . Consequently, 
for any given m one can construct [m/2 -|- IJ inequivalent spin Calogero models of Ujv 
type with PSRO, where [-J denotes the integer part. Among these models, only those 
with mi = m 2 (for even m) or mi = m 2 -I- 1 (for odd m) are unitarily equivalent to the 
su(m) spin Calogero model of Hat type with standard spin reversal operators introduced 
in Ref. [39] . 

As is t he case with the latter model, the configuration space C of the Hamilto¬ 
nian (2.15) can be taken as one of the maximal open subsets of on which the linear 


functionals xt ± Xj have constant signs. We shall again take C as the principal Weyl 
chamber of the root system, namely 


C = {x e : |a;i| < a;2 < • • • < xn } 


(2.17) 


Note that this configuration space contains its BCn counterpart (2.13) as a subset. As 


before, the Hamiltonian (2.15) is defined on a suitable dense subspace of the Hilbert 
space L‘^{C) 0 S. 

We shall next explain in detail how to construct the Uiv-type PF chain with PSRO 


associated to the spin dynamical model (2.151 by means of Polychronakos’s freezing trick. 
To begin with, note that the Hamiltonian g^n be decomposed as 




where 


i¥=3 


.K-)^ (4)^ 

is the Hamiltonian of the scalar Dn Calogero model and 




1 + Pjj 




1 + P 


{mi,m2) 


( 4 )" 


(2.18) 


(2.19) 


( 2 . 20 ) 
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is a spin-dependent multiplication operator. On the other hand, in the strong coupling 
limit a —>■ 00 the coefficient of the dominant term (of order a^) in the Hamiltonian (2.15) 
is given by 

1 1 

+ 


C/(x) = 








^2 


( 2 . 21 ) 


Hence as a —> oo the particles concentrate at the coordinates of the unique minimum 
^ of the potential I7(x) in the configuration space C [IS], and t he coo rdinate degrees of 
freedom of decouple from the internal ones. By Eq. (2.18), in this limit the 

eigenvalues of approximately given by 


E,, ~ Ef + a£j, 

where Ef^ and Sj are two arbitrary eigenvalues of and 


J^{mi,m2) — ,m2) _ 'y 




1 




( 2 . 22 ) 


(2.23) 


We shall take Eq. (2.23) as the precise definition of the Hamiltonian of the D^-type 


PF chain with PSRO. In fact, using Eqs. (2.10) and (2.11), it is easy to show that the 
Hamiltonians and are related by 


rp^^{mx,m2)rp _ -^(m2,mi) 


(2.24) 


Thus, we may assume without loss of generality that mi ^ m 2 , so that there are again 
[m/2 -I- IJ inequivalent PF chains of type with PSRO. Since the sites of these chains 
depend only on the scalar potential (2.21), the above models reduce to the su(m) PF 
chain oi Dm type with standard spin reversal operators |39j when mi = m 2 (for even 
m) or mi = m 2 -I- 1 (for odd m). See, e.g., Fig.j^for a comparison of the spectra of the 


Dm chain with PSRO (2.23) with mi = 3, m 2 = 1 and the su(4) DM-type PF chain with 


standard time-reversal operators (corresponding to mi = m 2 = 2) for = 10 spins. 

A brief remark on the relation between the Dm and BCm spin chains with PSRO in 
Eqs. (2.23) and (2.14) is now in order. As shown in |53|, the lattice sites of the former 
chain are given by ^1 = 0 and = -y/2i/i_i (2 ^ i < N), where yt > 0 denotes the 
A:-th root of the generalized Laguerre polynomial A)y-i- From the well-known identity 
NLjj^{y) = —yL]^_-^{y) and the previous characterization of the sites Q of the BCm 
chain ( 2.14| ), it immediately follows that ^ = lim^_>o C- Although one may naively 
think that the Hamiltonian jg giniply the /? —>■ 0 limit of its BCm counterpart 


n 


(rni,m 2 )^ this is Certainly not the case. The point is that, although the roots C,i with 
2 ^ i ^ N tend to finite nonzero limits when /3 —>■ 0, the first root Ci tends to 0 in this 
limit. As a consequence, the z = 1 term of the last sum in Eq. (2.14) need not vanish as 
/? —>■ 0, and in fact it can be shown [32] that 


lim -^ = Cr- 
/3^o Cf 


N 

2 


(2.25) 


Letting /3 —>■ 0 in Eq. (2.14) and using the latter identity we immediately obtain 
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Figure 1: Degeneracy d (in units of 10^) versus energy £ of the Djv chain l |2.23[ ) with mi = 3, m 2 = 1 
(blue), compared to the su(4) Djv-type PF chain with standard time-reversal operators (red), for N = 10 
spins. 

Thus, the /3 —)■ 0 limit of the Hamiltonian differs from its I?jv counterpart 

-^(mi,m 2 ) surface term or impurity interaction iV(l — /2. It is easy to 

see that this term vanishes only for e = 1, mi = m, m 2 = 0 (or, equivalently, e = — 1, 
mi = 0, m 2 = m). For any other choice of mi and m 2 , this surface term is nonzero and 
does not commute with the Hamiltonian Thus, except in the previously noted 

special cases, the spectrum of cannot be obtained from that of by 

taking the /3 —>■ 0 limit. This fact is illustrated in Fig. which shows that the spectra 
of these chains with mi = 3, m 2 = 1 and = 10 spins are clearly different. 



Figure 2: Degeneracy d (in units of 10^) versus energy £ of the Djv chain l |2.23[ ) with mi = 3, m 2 = 1 
and = 10 spins (blue), compared to its BCn counterpart in Eq. l |2.14[ ) (green; recall that the spectrum 
of the latter chain does not depend on /3) 


8 






















































































3. Spectrum and partition function 


In this section, we shall compute in closed form the spectrum and partition function 


of the spin Calogero model of Dm type with PSRO in Eq. (2.151. This will enable us to 
compute the partition function Qf U^-type PF chain with PSRO (2.231 by 


a standard freezing trick argument. Indeed, from Eq. (2.221 it is straightforward to derive 
the following exact formula for in terms of the partition functions and 


Z of the spin dynamical model (2.15) and of its scalar counterpart (2.19): 


Z('"i’™=^)(T) = lim 

Since Z has already been computed in Ref. [I 
evaluating once jg known. 


Z{aT) ' 


(3.1) 


, Eq. (3.1) provides an effective way of 


The key idea for deriving the spectrum of the spin Hamiltonian (2.15) is to ob¬ 


serve that it can be obtained by applying a suitable projection to a simpler differential- 
difference operator H' acting on scalar functions. The spectrum of H' can be readily 
computed by constructing a (non-orthogonal) basis of its Hilbert space on which this 
operator acts triangularly. The spectrum of is then easily determined by pro¬ 

jecting onto the Hilbert space of the latter operator. 

More precisely, the auxiliary operator H' is given by [39] 


I * 




(3.2) 


where Kij and Ki are coordinate exchange and sign reversing operators, defined by 

f (^1 5 ■ • ■ 5 ; • ■ • ) ^ ^ ^n '} f J 5 ■ ■ • 7 ^ ^ ^n') ^ 

K^f{xi,...,Xi,...,XN) = f{xi,...,-Xi,...,XN), 

and Kij = KiKjKij. The domain of the operator H' is of course a suitable dense subset 
of the Hilbert space The operator H' can be expressed in terms of the HAr-type 

rational Dunkl operators [15] 


J~ = 


dxi 


r,j¥=^ 


^^(1 — Kij) H—^7(1 — Kij) 
Xij Xjj 


(3.3) 


as m 


where 


H' = p(x) 




P(x) , 


(3.4) 


"2 


p(x) = e 4’- 

i<j 

is the ground state of the scalar Calogero model of D^-type and 

Eo = Na{a{N-l) + i) 
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(3.5) 


















is its ground-state energy. A basis of this Hilbert space on which H' acts triangularly is 
provided by the functions 


0„(x) = p(x)]^x"S n= (ni,...,nAr), (3.6) 

i 

where the rii’s are arbitrary non-negative integers. Indeed, since J~ lowers the degree 
|n| = ni + ■ ■ ■ -\- riN of any monomial from Eqs. (3.3) and (3.4) it immediately 

follows that 

H' (l)n{:x.) = E'^(l)n{^) + ^ Cmn</>m(x), (3.7) 

|m|<|n| 

where the coefficients Cmn are real constants and 

£;:, = a|n|+£;o. (3.8) 

As the diagonal elements of any upper triangular operator coincide with its eigenvalues, 
the spectrum of H' is given by Eq. (3.8). 

The spectrum of the spin Hamiltonian derived from that of H' by 

noting that these Hamiltonians are formally related by 






(3.9) 


In order to take advantage of this observation, we introduce the operator pro¬ 
jecting the Hilbert space 0 S onto states that are antisymmetric under particle 

permutations and symmetric under the action of ^ ^ j 

In other words, the projector is determined by 




where 


so that 


'^ij — ; 




(3.10) 


(3.11) 


_ _p. . ^('mi,m2) ^.^.^(™l."i2) _ p('m-i,m2) p(mi,m2) j^(mi,m2) ^2) 

We shall now outline the construction of the projector in terms of the analo¬ 
gous projectors for the BCN-type spin Calogero model with PSRO (2.14) with 

chirality e = ±1 (cf. [I]). To this end, recall that projects from the Hilbert 

space L^(M'^) 0 S onto spin wavefunctions antisymmetric under particle permutations 
and with parity ±1 under j g^^ 


.(mi,m2) _ a("*i,™2) 

oiVo _L — ii-i: 


The projector can then be expressed as 

N 


^{m i,m 2 ) j^{m 1,1712 


B.± 


= ±A 


(mi,m2) 
B,± > 


. (mi,m2) _ _^ 

B,± - 2N]^\ 


N N\ 

n(i±,("..-))^^„j,,, 


(3.13) 


(3.14) 




1^1 
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where Vi denotes an element of the realization of the permutation group generated by 
the operators tt^- and £/ is the signature of Vi- From Eqs. (3.101 and (3.13) we conclude 
that 

(3.15) 


A (mi,m2) _ A (mi,m2) A (mi,m2) 


Indeed, the right-hand side of the latter equation is clearly a projector, since 


^(mi,m2)^(mi,m2) _ ^(mi,m2) ^(mi,m2) _ q 








and it satisfies (3.10) on account of (3.13). Thus the space 




decomposes as the direct sum 


V = Vb,+ © Eb.-, 'FB.i = 


(3.16) 


We have already mentioned that, due to the impenetrable nature of the singularities of 
the Hamiltonian j^g Hilbert space can be taken as the space L^{C) ©5 of spin 

wavefunctions square integrable on the open set C in Eq. (2.17). On the other hand, 
any point in not lying on the singular subset Xiztxj=Q, can be 

mapped in a unique way to a point in C by a suitable element of the Djsi Weyl group, 
which is generated by coordinate permutations and sign reversals of an even number of 
coordinates [15]. Using this fact, it can be shown that L^{C) © 5 is actually isomorphic 
to the space V, and is equivalent to its natural extension to the latter space 

which (with a slight abuse of notation) we shall also denote by With this 

identification, in view of Eq. (3.12) we can write 


jj(mi,m2) _ J^(mi,m2)^(mi,m2) _ jjf^{mi,m2) 


(3.17) 


where H' acts trivially (as the identity) on S. 

We shall now explain how the spectrum of i 7 (™i’™ 2 ) derived from that of H' 


using the previous equation. To this end, note that by Eq. (3.16) the Hilbert space V is 


the closure of the linear subspace spanned by the spin wavefunctions 


^^,,(x) = A^7’"*=)(<(.„(x)|s)), 


e = ±, 


(3.18) 


where |s) = |si,..., sn) is an arbitrary element of the canonical spin basis. In fact, the 
wavefunctions (3.18) with^xed e span a subspace whose closure is the Hilbert space Vb,£- 
Clearly, the functions (3.18) are not linearly independent. Indeed, using Eq. (3.13) it is 


easy to show that these functions satisfy the relations 

^n.sW = > V'n.sW = e(-l)”‘ + ^(®* 


W: 


(3.19) 


where n' and s' are respectively obtained from n and s by permuting any two of their 
components (the same for both). Due to these identities, the sets {V'is(^)} {V'ir.s(^)} 

are both linearly independent provided that the following three conditions are imposed 
on the quantum numbers n and s: 
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i) To avoid overcounting, and for later convenience, we shall order the components of 
n as follows: 


n = (no, no) = ( 2pi,..., 2pi, ..., 2p ,,..., 2ps, 

h h 

2<?i + 1) • ■ •) 2gi + 1, ..., 2qt + 1,..., 2qt + 1) j 

where 0 ^ s,t ^ N, pi > p 2 > ■ ■ ■ > Ps ^ 0 and qi > q 2 > ■ ■ ■ > qt ^ 0. 


ii) By the second equation in (3.191, the allowed values of Si corresponding to each nt 
are given by 


Si S 


{1,2,..., mi} , for even rii , 

(mi + 1, mi + 2,..., mi + m 2 } , for odd rii , 


for the set by 


Si G 


for odd rii 


(1,2,...,mi}, 

(mi + 1, mi + 2,..., mi + m 2 } , for even rii ■ 


for the set {'0n.s(x)}. 

iii) If rii = nj and i < j we shall take Si > Sj, again to avoid overcounting. 

If the above conditions are satisfied, each of the sets {'0n,s(^)} = ±) is a non-orthogonal 

basis of the corresponding subspace Ve.e, and the union of these sets provides a non- 


orthogonal basis of the whole Hilbert space V by Eq. (3.16). We shall next show that 
leaves invariant each of the subspaces Vb^e, and that it acts triangularly on the 
corresponding basis s(x)} provided that we (partially) order it by the total degree |n|. 

Indeed, using Eqs. (3.15) and (3.17), and taking into account that [H', = 0 we 

obtain 

= Ag{-™=)((H'<).„(x))|s)). (3.20) 


From this equation and Eqs. (3.7) and (3.18) it readily follows that 




^ ^ C’mn'^m.s' (x) > 

|m|<|n| 


(3.21) 


where the Cmn’s are real constants and s' is a permutation of s s uch that (m, s') satisfies 
conditions i)-iii) above; see [51] for more details. By Eq. (3.21), the action of 


on the whole Hilbert space V = Vb,+ © ke,- is the direct sum of two upper triangular 
actions on each of the subspaces Vb.±- Consequently, the eigenvalues of this operator are 
given by 

K,s = K = a\n\+Eo. (3.22) 


where e = ± and (n,s) satisfies conditions i)-iii) above. Since the RHS of Eq. (3.22) 


does not depend on e and s, the eigenvalue associated with the quantum number n has 
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an intrinsic degeneracy coming from the two possible chiralities and the spin 

degrees of freedom. This intrinsic degeneracy is in fact the sum 


“n “ “n.+ "T “n,- i 


(3.23) 


where is the number of spin states satisfying conditions i)-iii) for the given n 

and e. Using these conditions we readily obtain [T] 


“n,+ 



,(mi,7712) _ ,(m2,mi) 
“n,- “ “n,+ 7 


(3.24) 


and therefore 


Ami,m 2 ) 

u-n 



(3.25) 


Thus the spectrum of the idjv-type spin Calogero model with PSRO (2.151 is given by the 
RHS of Eq. (3.22), where each level possesses an intrinsic degeneracy given by Eq. (3.25). 
Of course, the actual degeneracy of an energy a|n| + Eq is the sum 




|n'| = |n| 


where the sum is over all multiindices n' satisfying condition i) above. 

It is worth mentioning at this point that the spectrum of the RCjv-type spin Calogero 


model with PSRO and chirality e in Eq. (2.5) is also given by the RHS of Eq. (3.22), 
with Eq replaced by [T] 

Eq,b = Eq + N. 

Moreover, the intrinsic degeneracy of the energy o|n|+ifo,B is given by It follows 

from Eq. (3.23) that the spin Hamiltonian jg constant) the direct 


sum of two HCjv-type spin Calogero models of opposite chiralities with PSRO. Using 


Eqs. (3.22) and (3.25), the canonical partition function of the ZlAr-type spin Calogero 


model with PSRO can be written as 


^ qEo/a E ^(mi,m2)^|n| 


g = e 


-iKkuT) 


(3.26) 


where the sum ranges over all multiindices n satisfying condition i) above. Similarly, the 
partition functions of the corresponding SCAr-type models (2.5) are given by 


7(mi,'m2) 

"B,± 


{aT) = g(^0,B)/a^^OM.-2)^|„| ^ ^(m^mi) 


(3.27) 


From Eq. (3.23) it then follows that 


,(mi,m2) 
'B 
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r(m2,r, 

"B,+ 


zrf’"‘^’(ar)]. 


(3.28) 












In order to apply the freezing trick formula (3.11, we need only recall the expression for 
the partition function Z of the scalar Calogero model of Dj^-type derived in Ref. |55] . 
namely 

q-^o/<^Z{aT) = (1 + q^) n(l “ + q^)Z^{aT ), (3.29) 

i 

where Zb denotes the partition of the scalar Calogero model of type. Dividing 

Eq. (3.28) by Eq. (3.29) and applying the analog of the freezing trick formula ( |3.1[ ) for 
the partition function pp chain of BCtq type (2.14) we finally obtain 


= (1 + 0”' + 


(3.30) 


where from now on we shall use the variable q = e place of T. The partition 

function turn be expressed in terms of the partition function (q) of 

the su(to) PF chain of type A with k spins with Hamiltonian 


n 


(m) 


= E 


1 + Pjj 


(3.31) 


where pi is the i-th zero of the Hermite polynomial of degree k. Indeed, it is shown in 
Ref. [1] that for TO 2 > 0 we have 


■'B,+ 


N 


( 9 ) = E' 


N-k 


k^O 


N 

-Nl 




(3.32) 


where the q-binomial coefficient [ J ^2 is defined as 


N 

k 


( 9 ") 


N 


{q^)k{q^)N-k ’ 




(3.33) 


i=l 


Combining Eqs. (3.30) and (3.32) we hnally arrive at the following expression for the 
partition function of the Djv-type PF chain with PSRO (2.23) in terms of its type A 
counterpart as 

N 

2(-i.-+( 5) = ^ Zi”;^)(g2) zi")^Lfc(9") (m 2 > 0), (3.34) 


k=0 


where fN,k{q) is given by 


fN,k{q) = 


qN-k qk 


1 + 9 


N 


J 


(3.35) 


The case m 2 = 0, for which = l and the Hamiltonian (2.23) reduces to the 

rational version of the (trigonometric) Simons-Altshuler chain 22 j deserves special at¬ 
tention. Indeed, in this case by Eq. (3.19) the components of the multiindex n are all 
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even (resp. odd) for the eigenfunctions (resp. g). As shown in Ref. [T], this entails 
that for m 2 = 0 Eq. (3.32) should be replaced by 




-B.H- ’(9) = 4"’(9"). 


(3.36) 


On the other hand, since = —1 we have 


.1,(0,'m 2 ) _ nj(mifi) 
'^B,+ ~ '^B,+ 


C2 — '^B.S- 


■N 


by Eqs. (A2)-(A5) of Ref. [35]. From (3.36) it then follows that 


(3.37) 


(3.38) 


and substituting into Eq. (3.30) we finally obtain 

Z(--°)(Q)=zi”)^)(g2). 

Note that, as shown in Ref. [I], the RHS of the latter equation also coincides with the 
partition function of the SCAr-type chain (2.14) with e = 1 and m 2 = 0. This was to 
be expected, as the latter model reduces to its Ojv counterpart (2.23) when m 2 = /3 = 0 
and its spectrum does not depend on /3. 

As is well known, several equivalent closed-form expressions for the partition function 
of the Afc_i-type PE chain (3.311 exist in the in the literature |3| [HI |44l|49|. For instance, 
Polychronakos |3] showed that this function is given by 

5 E fci(fci-i) 


( 9 ) = 9 

fciH- \-km. — k 


[/ci, . . . , A^m] 


<1 5 


(3.39) 


where the q-multinomial coefficient [fci, 

[ki, 


^m\q 


is defined by 




(9) 


fclH-|-fcT7 

m 

n {q)ki 


Another well-known expression for the partition function was derived in Ref. m, 
namely 


E -^3 


k—r 


11(1 


(3.40) 


feVk 1=1 

Here Vk represents the set of all ordered partitions f = {/i, / 2 , ■ • ■, fr} of the integer fc, 
dm(f) = nl=i (y)) fi til® partial sums of f, and the complementary 

partial sums are defined as {E(, ,..., = {1, 2,..., fc} \ {Ei, J 2 , ■ ■ • i-^r}- A 

related expression for the partition function of the chain (3.31) can be obtained by 
exploiting its connection with a one-dimensional classical vertex model consisting of 
A: -I-1 vertices connected by k intermediate bonds [44]. Any possible state for this vertex 
model can be represented by a path configuration given by 


S — {^l, S2, . . . , Sk} • 
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(3.41) 










where Si € {1,2, •• • ,m} denotes the spin state of the i-th bond. The energy function 
associated with this spin path configuration s is defined as 

fc-i 

= (3.42) 

J=i 

where 0 is Heaviside’s step function, defined as 


0(x) = 


0 , if a; < 0 , 
1, if a; ^ 0 . 


(3.43) 


Using the Yangian quantum group symmetry of the model (3.311, it can be shown that its 


partition funct ion coincides with that of the one-dimensional vertex model with energy 
function (3.42) (cf. 


]). Thus Z^f((q) can be expressed as 






(3.44) 


where the sum has runs over all possible spin path configurations. In particular, from 


Eq. (3.44) it follows that 


4 %) = 9 ^ 


_ ^ k(k — l) 


Thus the partition function (3.34) with m 2 = 1 reduces to 


N 




(3.45) 


k=0 


It is obvious from any of the expressions (3.39), (3.40) or (3.44) that the partition func¬ 


tion Z^J(q) is a polynomial in q. In particular, from Eq. (3.38) it follows that 

is an even polynomial in q, and its energies are therefore even nonnegative integers. By 


Eq. (3.34), to show that the partition function of the H^y-type PF chain with PSRO 
is a polynomial in q when m 2 > 0 it suffices to prove that the coefficients fN,ki<l) in 
Eq. (3.35) depend p olyno mially on q. Although it is well known that the g-binomial 
coefficient [^] ^ in (3.33) is indeed an eyen polynomial in q of degree 2k(N — k) |S0], 


it is not clear whether fN,k(<i) is also a polynomial. In fact, we haye yerified that this 
is the case for a wide range of yalues of N and all k ^ N. We conjecture that this 


is true in general, so that when m 2 > 0 the energies of the spin chain (2.23) are also 


nonnegatiye integers. Note that the latter fact also follows from the freezing trick for¬ 


mula (2.22), Eq. (3.22) for the spectrum of the spin dynamical model (2.15) and the 


analogous formula for the scalar Hjy-type Calogero model in Ref. 


4. Statistical properties of the spectrum 

A characteristic property of all spin chains of Haldane-Shastry type is the fact that 
their leyel density approaches a Gaussian distribution as the number of spins tends to 
infinity. This property has been rigorously proyed for the chains of type and their 

16 

















related one-dimensional vertex models and has been numerically checked for 

the BCj\! and Dj^ type chains with standard spin reversal operators I37H3I1I1I]. 
More recently, it has been established that the level density of the BCN-type PF chain 
with PSRO shows a similar behavior [T]. It is therefore of interest to ascertain whether 
the level density of the Z^Ar-type spin chain with PSRO in Eq. ( |2.23[ ) becomes normally 
distributed as the number of spins tends to infinity. In fact, Figs. [^ andclearly suggest 
that this is actually the case. We shall restrict ourselves in the rest of this section to the 


case 1712 > 0, since for m 2 = 0 the spectrum of the chain (2.23) is twice that of an su(mi) 
PF chain of v4jv_i type (with the same degeneracies) on account of Eq. (3.38). 

The spectrum of the spin chain (2.23) can be determined for any fixed N by evaluating 


its partition function (3.34) with the help of, e.g., Mathematica. It turns out that the 


most efficient way to compute the partition function appearing in the latter equation 
is using the recursion relation 


j.(m) 

-'A.fe 


min(m,A;) 


(9) = 




Z-1 


^k-l 


nu 


^k—i 


)-4”L(9) 


(4.1) 


with the initial condition (q) = I (see Appendix A). In this way it is possible to 
evaluate the partition function qj^ ^ standard desktop computer for relatively 

high values of N (of the order of 50) and, say, mi -I-TO 2 < 4. Our computations show that 
the energy levels of the D^-type spin chain with PSRO are always a set of consecutive 
integers. This result is consistent with the fact that the spectrum of all previously studied 
rational spin chains of HS type is a set of consecutive integers [3l |38l |39] , including the 
rational spin chain of BC^ type with PSRO introduced in Ref. [T]. For this reason, in 
order to test the Gaussian character of the level density of the chain (2.23) as N —>■ 00 
one can compare directly its normalized level density 


fi£)=m-^^d,S{£-£,) 


m = mi + m 2 , 


(4.2) 


i=l 


where £i < ■ ■ ■ < £l are the distinct energy levels and di is the degeneracy of £i, with 
the Gaussian distribution 

1 

g{£) = ^^e-^ (4.3) 

with parameters /i and tr given by the mean and standard deviation of the spectrum, re¬ 


spectively. More precisely, the level density of the chain (2.23) is asymptotically Gaussian 
provided that 

di 




m 


A^» I. 


In order to check the validity of the latter equation for any given mi, m 2 and N we 
need to compute the corresponding values of /i and a. We shall next show that, as is 
the case with other spin chains of HS type, these parameters can be easily evaluated in 
closed form from their definition 




= m-^ 


cr^ = m“^tr 






(4.4) 
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The traces appearing in (4.4) can be computed in essentially the same way as for the 
BCjv-type PF chain with PSRO (2.14), using the traces of the spin operators Pij, 


L LJ 


given in Ref. [T]. Proceeding in this way we obtain 




(4.5) 




a^ = 


^ ■*" ^2 (4-6) 


where t = mi — m 2 and 


^ij ~ ij) 


hij — (Ci + Ci) 


The sums in Eqs. (4.5)-(4.6) can be evaluated by taking the /3 —>■ 0 limit of the corre¬ 
sponding formulas in Appendix A of Ref. [SB]- We thus obtain 


= + )«(«-!), 


1 


36 


1 


= hF (1 - f:2 ) ^(^ - 1)(4^ + 1) + - 1). 


(4.7) 

(4.8) 



Figure 3: Left: level density of the chain l |2.23| | with mi = 3, m 2 = 1 and N = 10 (blue dots) compared 
to the Gaussian distribution ( |4.3[ l (continuous red line). Right: analogous plot for W = 20 spins. 


We have checked that the normalized level density of the spin chain ( |2.23[ ) is indeed 
in excellent agreement with the Gaussian distribution (4.3) for different values of mi, 
m 2 , and even moderately large values of iV > 15. As an example, in Fig. we compare 
the normalized level density of the chain (2.23) with mi = 3, m 2 = 1, for iV = 10 
and A = 20 spins, respectively, with the corresponding Gaussian distribution (4.3). 
It is apparent from these plots that the fit, already quite good for N = 10, improves 
significantly for N = 20. This is confirmed by computing the RMSE errors for both fits, 
which are respectively equal to 3.66 x 10“^ and 2.18 x 10“^. Eor comparison purposes, 
we note that this error decreases to 1.11 x 10 ^ for A = 50 spins. 

Another interesting property of the spectrum of the chain (2.23) is connected to the 
distribution of the spacings between consecutive levels of the unfolded spectrum |53], 
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which in this case is given by 


Si = iVi+i - 


z = 1,... ,L - 1, 


where A = {rjL - m)/{L - 1), ??i = and 


rjiS) = f g{£')d£' 

J —OO 


1 

2 


1 + erf 



According to a celebrated conjecture due to Berry and Tabor [^, the distribution of these 
spacings for a “generic” quantum integrable system should be Poissonian, i.e., p{s) = e“®. 
On the other hand, a fundamental conjecture in quantum chaos due to Bohigas, Giannoni 
and Schmit [55) posits that the spacings distribution for a fully chaotic quantum system 
invariant under time reversal should follow Wigner’s law 

p{s) = (7rs/2)exp(—7rs^/4), 


characteristic of the Gaussian orthogonal ensemble in random matrix theory |56j . In fact, 
it has been shown that the spacings distribution of a large class of integrable spin chains 
of Haldane-Shastry type follows neither Poisson’s nor Wigner’s law m m m EH EH. 
More precisely, it is shown in Refs. mm [58] that the cumulative spacings density 
P{s) = /p p(s')ds' of a quantum system with equispaced energy levels and asymptotically 
Gaussian level density follows the “square root of a logarithm law” 


P{s) 1 - 


/TTSn 




£l-£i 

\f^ a 


(4.9) 


provided that a few mild technical conditions are satisfied. We have just shown that 


the energy levels of the rational Dm chain with PSRO (2.23) are equispaced and its 


level density is asymptotically Gaussian, and it can be easily checked using the formulas 
for £i and £l below that the technical assumptions in Ref. [58] are satisfied. Thus the 


spacings distribution of this chain is again approximately given by Eq. (4.9). It should 


be noted that for a more precise test of the validity of the Berry-Tabor conjecture one 
should restrict oneself to eigenspaces with well-defined quantum numbers corresponding 
to the main symmetries of the model. On the other hand, the fact that the spacings 
distribution of the whole spectrum is not Poissonian suggests that the Berry-Tabor 
conjecture does not hold in these eigenspaces, since the superposition of even a small 
number of Poissonian distributions is also Poissonian [59] . 

One of the characteristic properties of both the original Haldane-Shastry and the 
Polychronakos-Frahm spin chains of A^v-i type is their invariance under the quantum 
group E(s1(to)). From the existence of such a large symmetry group one should expect 
that the spectrum of these chain exhibits a high degree of degeneracy. In fact, it is 
shown in Ref. (50] that the spectrum of these models is far more degenerate than that 
of a generic Yangian-invariant system, due to their equivalence to a vertex model of 


the form (3.42) with a very simple dispersion relation. Indeed, as shown in the latter 
reference, the number of distinct levels of a generic Y(sl(TO))-invariant spin system 
with a large number of sites N behaves as where 1 < Am < 2 is the highest real root 
of the polynomial A’" — A"*“^ — • • • — 1. In contrast, grows as a polynomial in N for 
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all spin chains of HS type associated with the ^Ar_i root system. For instance, in the 
case of the type Ajv-i PF chain this polynomial is simply given by — £i + f, since its 
spectrum is a set of consecutive integers. From the explicit expressions for the maximum 
and minimum energies of this model in Ref. m we easily obtain 


- 1 


m ) 


l{m — 1) 
2m 


+ 1 


(PF chain), 


(4.10) 


where I = iVmodm. The situation is far less clear for spin chains of HS type associated 
to other root systems, with either standard or polarized spin reversal operators. On the 
one hand, the presence of these spin reversal operators breaks su(to) invariance, so that 
it is not obvious whether these models are invariant under a suitable quantum group, let 
alone T(s1(to)). On the other hand, it has been observed that the spectrum of some of 
these chains is also highly degenerate, which seems to indicate the presence of a large 
symmetry group. 


In the particular case of the D^-type chain with PSRO in Eq. (2.23), the number of 


distinct energy levels can again be exactly computed under the assumption (which we 
have numerically checked) that the spectrum consists of consecutive integers. Indeed, it 
suffices to evaluate the maximum and minimum energies ^ and in terms 

of which the number Qf distinct energy levels is given by 

(mi,m 2 ) c(mi,m2) . -i 

^ ^mav • 


In the first place, the maximum energy can be easily computed by taking into account 
that Pij and self-adjoint operators whose square is the identity, so that 

their eigenvalues are ±1. Moreover, it is clear that a state of the form |s, s,..., s) is a 


simultaneous eigenvector of all the operators Pij and P-^ 


(ml,m 2 ) 


the maximum energy of the chain (|2.23 1 is given by 


with eigenvalue 1. Hence 




-2 




(4.11) 


where the sum was evaluated in j39j . On the other hand, by Eq. (3.301 the minimum 
energy is given by 

( 4 . 12 ) 


Cin - mm I 


where minimum energy of the BC^-type chain (2.14) with e = -1-1. The 

latter energy was computed in Ref. [T], with the result 


4“+’m2) = {N- 1){N + I - mi)/m + {I - mi)9{l - mi) 


(4.13) 


where I = Nmodm and 9 is Heaviside’s function (cf. Eq. (3.43)). Using the above 
relation it is straightforward to check that if mi ^ m 2 we have 4+ ^ 4”+’™^t and 

therefore 

mi ^ m 2 . (4-14) 


p(mi,m2) _ 

^min ~ ^B,+ ’ 


From Eqs. (4.11) and (4.14), and the assumption that the energy levels are equispaced, 


we finally obtain the following closed formula for the number of distinct energy levels of 
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the D]\f chain (2.23): 


m J m 


m 2 ,, , 1(1 - mi) 


— (I — mi)9{l — mi) + 1 . 


(4.15) 


Thus, it is apparent that jg quadratic polynomial in TV, as is the case with the 

PF chain of Ajq-i type (cf. Eq. ( |4.10 )). In particular, the spectrum of the chain (2.23) 
exhibits a very high degeneracy, much larger than that of a generic Yangian-invariant 
su(m) spin model; see, e.g., Fig. 



N 


Figure 4: Lo garithmic plot of the number of distinct energy levels of the rational D^v-type chain with 
PSRO (|2.23|l with nii = 3, m 2 = 1 (red dots) and of a generic su(4) Yangian spin model (red dots) for 
10 ^ 


The low number gf distinct energy levels of the model (2.23) entails an 

extremely high average degeneracy i^{mi,m 2 )^ which in turn suggests the 

existence of a large symmetry group. More precisely, it was shown in Ref |60j that the 
polynomial growth of the number of distinct energy levels of the spin chains of HS type 
associated to the A^-i root system is ultimately due to the equivalence of these chains to 
a Yangian-invariant vertex model of the form (3.42) with a suitable dispersion relation. 
This observation makes it reasonable to conjecture that the D]\f-type spin chain with 
PSRO (2.23) is also invariant under a suitable Yangian group, and that its spectrum 
coincides with that of a vertex model analogous to (3.42) with an appropriate energy 
function. 


5. The ferromagnetic models 

We shall consider in this section the ferromagnetic counterparts of the D]\[-type spin 
Calogero model with PSRO (2.15) and its corresponding spin chain (2.23), with Hamil¬ 
tonians respectively given by 




{mi,m 2 ) 


L 


(X ^ij 


_ p(mi,m2) 


(4)^ 


2 


(5.1) 


^(mi,m2) 

rip 


E 


1 - Fp 
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1 _ p(”^l’'^2) 


(Ci + 


(5.2) 


and 


































The spectrum of the ferromagnetic spin Calogero model (5.11 can be studied in a similar 


way as its antiferromagnetic counterpart, following the procedure described in Section]^ 
To begin with, we note that the Hamiltonian (5.1) and the auxiliary operator ( |3.2[ ) are 
related by 






(5.3) 


Hence, the operator in Sectionshould be replaced by the projector 

onto states symmetric under simultaneous exchange of the particles’ spatial and spin 
coordinates, and with parity +1 under the product of an even number of operators 
^(mi,m 2 ) (3.11)). The new projection operator is the sum 


A (mi,m2) _ A (mi.m2) I A (mi,m2) 

of the symmetric analogs of the HCAi-type projectors in Section]^ determined by 


.(mi,m2) _ .(mi,m2) 

l-Zl-'D ^ -L. - 


(mi,m2) A (mi,m2 
f..- 


_A (mi,m 2 ) 


(5.4) 


’’’b^''B,s,± ~ ^'-B,s,± ’ '4 ^‘-B,s,± 

As explained in Sectionj^for the antiferromagnetic case, the operator jg equiv 

alent to its natural extension to the Hilbert space 


H = Hb,s,+ © Hb,s,-, h^B,s,± = ® S)- 


(5.5) 


A set of (non-orthogonal) vectors whose linear span is dense in each of the Hilbert 
spaces Vb,s,± can be constructed in much the same way as in the antiferromagnetic 
case, replacing by in (3.18). Due to the symmetry of under 


permutations, in order to obtain a basis of these Hilbert spaces we must replace condition 
iii) in Section]^ by 

iii') Si ^ Sj if rzi = nj and i < j ■ 


As a result, the spectrum of the ferromagnetic model (5.1) is still given by Eq. (3.22) 
but the corresponding degeneracy factor in (3.25) should be replaced by 


j(mi,m2) 

-‘F.n 


n 


-\- ki — 1 
ki 


n 

i=i 


fm2 + lj - I 

V 


n ( W 2 + fci — 1 

V h 

l—L ^ 


t 

n 

i=i 


TOl + Ij — 1 

u 


(5.6) 


Using this formula for the degeneracy factor and proceeding as in Section we find 
that the partition function of the D^v-type ferromagnetic spin chain (5.2) is given by the 
following analog of Eq. (3.30[): 




-1 


,.(mi,m2) , 


-'B,F,+ (?) ^B,F,+ 




(5.7) 


where denotes the partition function of the ferromagnetic counterpart of the 


rational HC^-type chain (2.14). Proceeding as in Ref. [T] one can readily prove the 
ferromagnetic version of Eq. (3.32), namely 


j,(mi,m2) 

■'B,F,+ 


N 


(9) = J2< 


N-k 


k =0 


7(™i) ('„2\ 2 
-'A,F,fcW ) ^A,F,N-k 

22 


-fe(9^) (m2>0). 


(5.8) 


























Here f. denotes the partition function of the ferromagnetic version of the su(m) 
PF chain of type Ajsi-i (3.31) with k spins, obtained replacing Pij by —Pij in the latter 


equation. Finally, from Eqs. (5.7) and (5.8) we immediately obtain the following explicit 


formula for the partition function of the ferromagnetic chain (5.2): 


N 


^ ^ (m2 > 0), (5.9) 




where /jv,/c(q) is again given by (3.35). For m 2 = 0, proceeding exactly as in Section]^ 
we obtain 

= (5.10) 

Several explicit expressions for the partition function Z^^^ ^ of the ferromagnetic PF 
chain of A]^_i type appearing in the previous formulas are known in the literature. The 
first of these expressions is the analog of Eq. ( |3.39| ) , namely 

?(™) 


'^A,F,fc('?) — [kl, ■ ■ ■ Tkm]q ■ 

fciH- \-km — k 


Alternatively, z''^^ ^ may be obtained from Eq. (3.40) replacing dm{i) by its ferromag¬ 


netic version dF,m(f) = 111=1 ^)- Finally, k i® ^1®° given by the RHS of 

Eq. (3.44) with 9{x) replaced by 1 — 6{x) in the definition (3.42) of From any 

of these explicit formulas for 2a,F, k{q)^™'\ it follows that this function is a polynomial 
in q. By Eqs. (|5.9|)-(5.10l the same is true for the partition function of the chain (2.23), 


provided that the coefficient fN,k{Q) is a polynomial in q. 

As is well known, the partition functions of the Ajv-i-type ferromagnetic and antifer¬ 
romagnetic PF spin chains satisfy a certain duality relation [3l [611162]. In fact, a similar 
relation also holds for PF chains associated with other root systems [H |38l [39]. In order 
to establish a duality relation between the partition functions of the ferromagnetic and 
antiferromagnetic spin chains of type with PSRO, it suffices to observe that their 


Hamiltonians (5.2) and (2.23) are related by 


n 


{mi,m2} 


n 


(mi,m2) 


= 2E = NiN - I) 


(5.11) 




(cf. Eq. (4.11)). This obviously implies that the eigenvalues of q^{mi,m 2 ) 


are also related by (5.11), so that their partition functions satisfy the duality relation 


2.^mi.m2)(^) ^ ^Ar(Ar-l)2-(mi,m2)(-^-l)_ 


(5.12) 


6. Conclusions and outlook 

We introduce the Dm spin Calogero model with PSRO and its associated spin chain 
of HS type, namely the Dm PF chain with PSRO. We solve the former model by finding 
a suitable (non-orthonormal) basis of its Hilbert space on which its Hamiltonian acts 
triangularly. From the spectrum of this model we are able to compute its partition 
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function in closed form, which yields the partition function of the spin chain via Poly- 
chronakos’s freezing trick. More precisely, we show that the latter partition function 
can be expressed in terms of the partition function of the type-A PF chain. Since the 
type-A partition function can be efficiently evaluated using a simple recursion formula 
that we also derive in this paper, we are able to exactly compute the spectrum of the 
Z? 7 V-type chain for relatively high values of N. In this way, we are able to study several 
global properties of the spectrum of the latter chain. In particular, we provide strong 
numerical evidence showing that its energy levels are a sequence of consecutive integers, 
and that its level density becomes normally distributed when the number of spins tends 
to infinity. From these facts we conclude that the spacings between consecutive levels of 
the unfolded spectrum follows a “square-root-of-a-logarithm” distribution, characteristic 
of most spin chains of HS type. We also determine the number of distinct energy levels 
of the spin chain, showing that it is a second-degree polynomial in N, as is the case with 
the PF chain of A^r-i type. For spin chains of HS type related to the Am-i root system, 
it is known EO] that the polynomial growth of the number of distinct levels is a conse¬ 
quence of the fact that these models are equivalent to a Yangian-invariant vertex model 
with linear energy function and polynomial dispersion relation. Our results strongly 
suggest that this is also the case for the present model, a conjecture which certainly 
deserves further study. In particular, the validity of this conjecture would also point out 
at the existence of a suitable Yangian symmetry for both the Z^Ar-type spin chain and 
the spin Calogero model with PSRO, as is the case with the rational and trigonometric 
Calogero-Sutherland models of A^r-i-type and their associated spin chains. 

The present work suggests some possible future developments. Among them, the 
most natural one would be to address the extension of our results to the Sutherland 
(both trigonometric and hyperbolic) models of BCn, Hat and Dn type and their related 
spin chains. From a more mathematical standpoint, the fact that the chain’s spectrum 


consists of integers leads us to conjecture that the function fN,ki<l) in Eq. (3.35) is a 
polynomial in q. Although this conjecture can be easily checked numerically, we have not 
been able to find an analytic proof thereof using the properties of q-binomial coefficients. 
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Appendix A. Recursion relation for the partition function of the PF chain 
of Afc_i type 


We shall provide in this Appendix a short derivation of the recursion relation (4.1) 
satisfied by the partition function of the su(to) PF chain of Ak-i type. The main 
idea behind the proof is to decompose the multiindex f S Pfe in Eq. (3.40) as 


f = (/i,...,/,_i,0 = (f,0, 

with 1^1^ min(m, k) and f € Vk-i- Setting s = r — 1 we have 

B'l — I = Bs—I -\- Bg—l -\- k — I . 
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(A.l) 





and therefore 


, J-U} = {K ■ ■■h-i-s] u {fc - / + 1,..., - 1} . (A.2) 

Substituting and ( |A.2[ ) into Eq. ( |3.40| ) we obtain 


4“’(9)= i: 

min(m,/c) 


1-1 






2=1 

Z-1 


^ E 711(1■ 4^/(9)’ 


fePfc-i 

(m) 


f=i V/* 


k—l — s 


H -S — 


^ n ( 1 - 4 ) 


i=l 


as claimed. As to the initial condition, from Eq. (3.40) with k = \ it easily follows that 
(?) “ From the recursion relation (4.1) with k = 1 we easily obtain (g) = 1. 
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